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ABSTRACT 

Simulations of large-scale structure formation predict the formation of sheet- 
and filamentary structures, which are often invoked as the origin of the Lya for- 
est. In their simplest description, these sheets and filaments require a differen- 
tial distribution of observed line-of-sight velocity widths (b) which will decrease 
as power-laws at velocities well above the observed peak in this distribution: 
for filaments, the differential distribution is dN/db oc 6^'^, while for sheets it is 
dN/db oc b^^. These functional dependences on b arise a priori due to the geom- 
etry of these absorbing structures - assuming random orientations relative to the 
line-of-sight - and are otherwise unrelated to the physical state in the absorbing 
structure. We find the the distribution at 6 > 35 in three previously published 
data sets to be steeper than dN/dB oc b~^ (99.99% confidence). This implies 
that evidence of the finite length of these kinds of absorbing structures is present 
in the 6-distribution data. 



Subject headings: cosmology: early universe — cosmology: intergalactic medium 
— galaxies: quasars: absorption lines 



Introduction 



The observational distribution of measured velocity widths (6) of Lya forest absorption 
lines may be parameterized as Gaussian, but with a significant non-Gaussian, high-velocity 
tail (> 55 km sec'^ cf. [Hu et al. 1995| ; |Lu et al. 1996| ; [Kim et al. T99^ ; [Kirkman fc] 
[lytler 1997| ; |Khare et al. 1997| for recent observational results). 
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Investigations based on spherical gaseous clouds bounded by external pressure ( Puncan 



|et al. 1991| ; [Petitjean et al. 1993| ) find that the 6-distribution can only be reproduced by a 
population with a strong correlation between central density, radii, and cloud mass, requiring 
a tight correlation between column density and b, which has not been strongly supported 
by the data QRauch et al. 19931 discusses selection effects which may have produced earlier 



reported correlations). Other researchers find that high implied temperatures of such clouds 



(assuming b = ^j2kT/mp) cannot be produced by structures in thermal equilibrium with an 
ionizing UV background, and may imply recent adiabatic collapse (Press fc Rybicki 1993|) . 

The ramifications of CDM structure formation models as the origin of the Lya forest 
have been extensively investigated. Recent workers, using hydrodynamic simulations of 
cosmological models (usually CDM or ACDM), describe the structures responsible for the 
low-density {Nh^IO^^) Lya forest to be filaments, pancakes, "cigar-like", mini-pancakes, or 
sheets (|Cen et al. 1994| ; [Zhang et al. 1995| ; [Hernquist et al. 1996| ; |Katz et al. T996| ; |Miralda-| 



[Escude et al. 1996| ; |Bi fc Davidsen 1997| ; |Cen fc Simcoe 1997| ; |Nath 1997t [Zhang et al.T997| ) 



Commonly, the 6-distribution is parameterized as a Gaussian; in addition, a non-Gaussian 
high-velocity tail is universally noted. The measured fe-values do not always correspond 
directly to a temperature, and most structure formation theories produce fiuctuations in the 
density of the inter-galactic medium (IGM) larger (in redshift space) than the thermal width 
( [Miralda-Escude et al. 1996[ ; [Hui et al. 19"97D . 



In this Letter, we describe how this high-velocity tail may be generated by the same 
structures responsible for the lower-velocity Gaussian peak. The low-velocity (Gaussian 
peak) values are caused by the line-of-sight passing nearly perpendicular to the long-axis 
of the sheet /filament-type structures, while the high- velocity tail is due to higher angles of 
incidence. The relative number of low vs. high 6-values are, in this scenario, determined by 
the geometry of the absorbing structures, and is otherwise unrelated to the physical state of 
the absorbing structure. 

In exploring this scenario, we keep in mind the fact that the fiuctuations which are 
produced in large-scale structure formation simulations are in a continuous density field, 
and ascribing a simplified geometry to them in the way we do (as a cylinder, or an infinite 
sheet) is an idealization of a structure which has already been observed in the simulations 
to be more complicated. 
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2. The Model: Line-of-Sight 6-Distribution 

In this model, the magnitude of b corresponds to the line-of-sight width of the absorbing 
structure in velocity space, which is determined both by its physical size in comoving space 
and by its peculiar velocity structure. 



2.1. Cylindrical Filaments 

In the simplest manifestation (refer to Fig. |I]), an observer's line-of-sight (LOS) crosses 
an absorber which has a minimum (velocity) width bo, and which has a long-axis. We assume 
that the geometry of the absorber is such that the cross section - which is the plane shared 
by the short-axis 6o and the LOS, parallel to the long axis - can be described as a rectangle 
of width bo and quasi-infinite length. This is a fair assumption to describe cross-sections of 
filaments and sheets, but not of spheres. Presently, we treat the long-axis as infinite, but 
the expected finite length would produce observationally important effects, which we discuss 
below. 

The LOS crosses the absorber at an angle 9 to the perpendicular of bo, producing a 
line-of-sight crossing distance (in velocity space) of b. The observed velocity-width b is 
geometrically related to the short axis bo: 



This is the relationship between the LOS velocity-width (6) of the observed absorption 
line and the velocity-width bo of the density enhancement (that is, filament or sheet). If we 
assume a random distribution in 6 (that is, a random orientation of absorbing structures 
relative to the LOS), because 6' = is the axis of symmetry for the cylindrical filament, the 
number of LOS per unit 6 is: 

sm{9) (2) 



and the number of LOS per unit velocity for a constant bo - is: 



dN, 



cyl 



db 



{b,bo 



dN de 
~dedb 



= b>bo 



(3) 



b < b 
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We assume a cylindrical filament, with a diameter L; bo then depends on an angle 
0, corresponding to an impact parameter relative to the center of the cylinder's circular 
cross-section: 



bo = Lcos{(j)) (4) 

where = when the LOS passes through the center of the cylinder's circular cross-section, 
and ~ 7r/2 when the LOS passes nearly tangent to the cylinder's circular cross-section. 
The resulting distribution of bo, assuming a random distribution of impact parameters, is: 



dN 



If we assume a random distribution of angles, the total 6-distribution, using Eqn. 
becomes 



(5) 



COS2(0) 

b> L 

arccos ij) b < L 



After making a change of variables of X = sin(0), converting the integrating variable 
from dcp to c?(sin(0)), and replacing j = Jb, the integral is in a tabulated form ( p wight 1961 



200.01 and 202.01). Once integrated, we find the distribution of velocity-widths b to be: 

)6 > 1 
)b < 1 
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This is the normalized 6-distribution function for an ensemble of cylindrical filaments of 
velocity-diameter L, oriented randomly with respect to the line-of-sight. This distribution 
function is plotted in Fig. ||a, super-imposed on combined datasets 1+3+5 (see Table |l|). 
The distribution is sharply peaked, with a broader wing at lower-velocities than at higher 
velocities. It is a very poor description of the data near the Gaussian peak, indicating that, 
at the very least, a distribution of L values must be imposed. 
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Expanding the logarithm in ]b in Eq. we find: 

plus higher order terms. Thus, the differential distribution of LOS velocities at values greater 
than the velocity diameter of the cylinder decreases as an inverse cube of the measured 
velocity. Because this decrease is a power-law, integrating over a distribution in L will not 
change it's power-law nature at 6 > Lmax (where Lmax is the largest contributing L value). 

The differential distribution can be integrated, producing the the normalized cumulative 
distribution: 

A, <^..n { \{W-{Jb-i))coth-im ]b>l 
'^-''^''•''•^'^I laO-(/^--J)tanlr'(«) ;,<1 

which is the fraction of h values expected to be larger than some value 6, for a known value 
of L. The median is /5 = 1.0 (that is, b = L). Approximately eighty per cent of the b values 
are between 0.5 and 2.0 of the median. At the high velocity end (/b ^ 1), the cumulative 
6-distribution decreases as /5~^, as expected. We show this cumulative distribution for an 
assumed L = 28 km sec~^, super- imposed on datasets 1+3-1-5 (see Table |lj), in Fig. |^. 
The power-nature of the cumulative distribution at b >28 km sec~^ is apparent, and is 
comparable to that seen in the dataset. (We comment on the possibility of a break in the 
data near 6 = 80 km sec~^ in Sec ^ where we statistically compare the observed distribution 
with a power-law distribution). 



2.2. Sheets 

The differential distribution for sheets may be found similarly to cylinders. The rela- 
tionship between b and bo (the short axis velocity) is as in Eq. |l], where again we refer to 
Fig. |l|. However, the axis of symmetry for the sheet is now perpendicular to the surface of 
the sheet, and the number of LOS per angle 9 changes accordingly (compare with Eq. |): 

^ = cosiO) (11) 

This is the important difference with the filament distribution. We assume a delta-function 
distribution of bo{=L - the width of the sheet). By chain rule, the differential 6-distribution 
is: 
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^(6,L)= Hf) (12) 
db ^ ' [0 b<L 

resulting in a cumulative distribution: 

iV,heet(>&;&,^) = j f (13) 

[1 b < L 

Thus, the differential distribution expected from sheets is oc and the cumulative 
distribution is oc b~^. The differential distribution is non-zero only at 6 >L - and thus is 
truncated at the sheet-width. 



3. Comparison with Observational Data 

The dependence of the b distribution dN/db oc b^^ or 6^^ may be present in existing 
data-sets. In this section, we perform a few comparisons between the data and power-law 
6-distributions. 

These comparisons may be hampered by contributions to line-measurements of line- 
blending (when 2 or more individual lines cannot be statistically decomposed, and thus the 
properties of only one, broader, blended line is measured), and also to the sensitivity of the 
detection measurement algorithms to lines of these velocity widths. We have assumed that 
the measurement algorithms are 100% efficient. 

We used data from three recently published studies of the Lya forest absorbers ( |Hu| 
et al. 1995| [Lu et al. 1996| ; [Kirkman fc 'lytler 1997| ). From these, we have used only lines with 



measured column densities Nh = 10^^'^ — 10^^, and 6 > 35 km sec ^. Using a KS-test ( [Press 



et al. 1995|) , we determine the probability that data from these observations are consistent 



with having been drawn from a differential distribution with a single-power-law, for a range of 
power-law values in 0.05 increments. We tested each data-set individually (datasets 2, 4, and 
6, see Table P, producing 99.5% confidence limits (where the dataset has a 0.5% probability 
of being drawn from a parent distribution with that power-law value). The results of these 
are shown in Table ||. The combination of all three datasets (dataset 7) places limits on a 
single-value-power-law parent population for the differential distribution with a power-law 
between [—4.05, —3.15]. The combined data-set is inconsistent with a single- value power-law 
dN/db oc at 99.99% confidence, and inconsistent with dN/db oc b'^ at (100-2 x 10-^9)% 
confidence. 
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In Fig. there is what looks hke an apparent break from a single power-law in the 
data near 6 = 80 km sec~^. In a combined dataset from all three observations, using values 
of 6 > 80 km sec~^, we find the 99.5% upper-limit on the power-law slope is —3.5 (the lower- 
limit is —11). This is consistent with the range found for 35 < 6 < 60 km sec~^( [-4.55,-2.05], 
99.5% confidence range, drawn from dataset 7; 346 lines). Thus, the break is not statistically 
significant. 

4. Discussion and Conclusions 

We have described a model in which we have assumed the absorbing structures of the 
Lya forest may be idealized as infinite sheets and cylinders in velocity space. In this model, 
the measured LOS velocity width b is due to the line-of-sight size in velocity space of the 
absorbing structure, determined both by its physical size in comoving space and by its 
peculiar velocity structure. Assuming a random orientation of these structures relative to 
the LOS, a high-velocity-tail to the distribution of LOS velocity widths of the Lya forest is 
produced. This tail is due to the chance occurrence of a LOS passing close to the long-axis, 
making the total traverse distance much longer than that of a LOS which passes nearly 
perpendicular to the long axis. The differential distribution due to an ensemble of infinite 
cylinders (sheets) shows a dependence A''cyi(> b) oc b~^ (iV'shoet(> b) oc b~^). 

The shape of the L-distribution (the characteristic size of the structure) will dominate 
the observational 6-distribution at low values [i.e. b ^ L). Thus, in this model the Gaussian 
which is observed in the velocity range b ^10-40 kmsec"^ is largely due to the shape of the 
L distribution, while the high-velocity tail is due to the line-of-sight penetrating absorbing 
structures close to their long-axis. 

The single power-law differential 6-distribution of b~^ for cylinders and for sheets 
should continue up to values of 6 ~ oo only if these structures are infinitely long. This is 
almost certainly not the case. At high values of b, the distribution should break from these 
power-laws, due to the finite size of the absorbing structures. 

A comparison with a combined dataset of 447 lines with b >35 km sec~^ places limits on 
the value of the power-law distribution, with dN/dB oc l)-['^-05,3.i5] (^99 5% confidence limits). 
A pure power-law oc 5~3.o jg excluded at a probability of 1 x 10~^ while oc is excluded at 
2 X 10~^^. Under this model, including values of b toward lower velocities (6 ) could 

bias the power-law toward more negative values; however, using a dataset consisting only of 
higher velocities {b > 80 km sec^^) the data are even more inconsistent with dN/db oc b^^. 
Therefore such low-6 contamination cannot be responsible for the inconsistency between the 
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observational data and the predicted power-law dependency. The effects of finite-length 
sheets and filaments may be responsible for the steepness of this power-law. Further work 
on the observational limits of such a break could provide a direct measurement of the size 
of these structures. Hydrodynamic simulations can play an important role in predicting the 
form and magnitude of this break, providing observers with the tools they need to find it. 

In conclusion, modeling the absorbing structures of the IGM responsible for the Lya 
forest as infinite cylindrical filaments and sheets in velocity space predicts power-law b- 
distributions - which predict more observed high-6 lines than are reported in the observations. 
This can be due to several causes, among them the finite-size of such structures, a decreasing 
detection efficiency with increasing b, or that the geometric idealization here presented does 
not accurately reflect the physical system. 
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fiir extraterrestrische physik. We are extremely grateful to Lam Hui, for very useful discus- 
sions and comments on drafts. This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, Cahfornia Institute 
of Technology, under contract with the National Aeronautics and Space Administration. 
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Fig. 1. — Schematic diagram of the absorption. The hne-of-sight crosses an absorber with 
a minimum size 60 and an infinite length, at an angle 9 with respect to the perpendicular, 
producing a total line-of-sight crossing distance of h km sec~^. 

Fig. 2. — Comparison of combined datasets 1+3+5 ( Table |l]) with theoretical 6-distribution 
due to a cylinder with a single radius L=28 km sec^^. (a) The differential distribution - the 
solid histogram is the data, and the broken line is the theoretical prediction, which is a very 
poor approximation to the data; this indicates that a distribution of cylinder radii is required, 
(b) Integral of the data (and theoretical model) shown in panel (a). It is shown in log-log, 
to make the power-law nature of both the theoretical distribution and the observed data 
clear. The apparent break near 6 = 80 km sec~^ in the data is not statistically significant 
(see Sec. |). 



Table 1. h Data-sets 



Dataset 


Nh range 


h range 


Number of 


dN/dh oc 6" 


# (ref) 


(cm~^) 


(km sec~^) 


h values 


(99.5% limits) 


1(A) 


IO12.5-14 




328 




2(A) 


IO12.5-14 


>35 


116 


[-4.1,-2.50] 


3(B) 


IO12.5-14 




790 




4(B) 


IO12.5-14 


>35 


252 


[-4.55,-3.05] 


5(C) 


IO12.5-14 




287 




6(C) 


IO12.5-14 


>35 


84 


[-4.70,-2.40] 


7(A+B+C) 


IO12.5-14 


>35 


447 


[-4.05,-3.15] 


8(A+B+C) 


IO12.5-14 


>80 


46 


< -3.7 



References. — (A) [Kirkman fc Tytler 1997| : QSO HS 1946+7658, (z=3.02) 
(B) |Hu et al. 19"95| : Q0014+813 (z=3.38), Q0302-003 (z=3.29) 0636+680 
(z=3.17), 0956+122 (z=3.30); (C) |Lu et al. 1996| : QOOOO-26 (z=4.1) 
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